Newman-Penrose quantities as valuable tools in 
astrophysical relativity 



Marco Bruni*, Andrea Nerozzi* and Frances White* 

* Institute of Cosmology and Gravitation, University of Portsmouth, Portsmouth POl 2EG, UK. 

Abstract. In this talk I will briefly outline work in progress in two different contexts in astrophysi- 
cal relativity, i.e. the study of rotating star spacetimes and the problem of reliably extracting gravi- 
tational wave templates in numerical relativity. In both cases the use of Weyl scalars and curvature 
invariants helps to clarify important issues. 



INTRODUCTION 

The Weyl scalars - the components of the Weyl tensor over a null tetrad - were known and 
used in relativity before the introduction of the Newman-Penrose (NP) formalism, but 
within the latter they acquired a new relevance. Here I will summarise the use of Weyl 
scalars as tools in two different contexts in astrophysical relativity. First I will briefly 
summarize recent work Jl|] aimed at assessing the validity of the Hartle-Thorne (HT) 
slow-rotation approximation for describing stationary axisymmetric rotating Neutron 



Stars (NSs), introducing work in progress [2] to extend the analysis. In this context the 
Weyl scalars are used to construct invariant measures of the deviation of the exterior 
spacetime from Petrov type D, in view of a possible development of a Teukolsky-like 
perturbative formalism for rotating NSs. 

In the second part I will outline how the Weyl scalars may be used in numerical 
relativity in order to construct a wave extraction formalism for simulations dealing with 
spacetimes that will settle to a perturbed black hole (BH) at late times. In this case the 
Teukolsky BH perturbation formalism [3] is in principle applicable, but it is difficult to 
extract a BH background spacetime, i.e. the gravitational mass and angular momentum, 
from a given simulation. Introducing the notion of a quasi-Kinnersly frame (also used 
in IH), in 0,0] a method was proposed that bypasses this difficulty, by not requiring a 
background, and allows direct wave extraction. I will present here work in progress Ja] 
where the method is directly applicable. 



THE SPACETIME OF ROTATING STARS 

Using a variety of models for NSs of different masses and equations of state, and 
comparing with full general relativistic numerical models, it was shown in |l|] that 
the HT approximation to the metric of a rotating relativistic star is very good for most 
astrophysical applications, even at the rotation rates of the fastest known milli-second 
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FIGURE 1. Contour plots of (1 - |5|)xl0 4 for decreasing values (from bottom left) 5,4,3,2, 1, for a 1.4 
M Q star, for a representative equation of state. Left panel: e = 0.39258; right panel: e = 0.54440. 

pulsar. For instance, the innermost stable circular orbit is predicted with an accuracy of 
~ 1%. It was also shown in [ 1] that although the spacetime of these stars is of Petrov type 
I (general; it would be type D for a spherical star), the deviation from type D is small, 
at least on the equatorial plane. The HT metric is obtained as a perturbative solution to 
second order in e = Q./CI*, where Q. is the star's angular velocity and CI* = (M/R 3 ) l / 2 is 
a "Keplerian" rotational scale. The deviation from type D is measured by 1 — \S\, where 
S = 27 J 2 /I 3 is an invariant curvature scalar, with S = 1 for type D. Fig.[l][2] shows that 
1 — \S\ is small also out of the equatorial plane. In fact, 1 — |5| decreases more rapidly to 
zero out of the equatorial plane, as one would have expected. Thus in this sense 1 — \S\ 
in the equatorial plane is a good upper limit for the deviation from type D. 



WAVE EXTRACTION 

In 0, 0] methods were introduced to identify, for a general numerical relativity imple- 
mentation, what was dubbed the Quasi-Kinnersley frame, an equivalence class of tetrads 
that reduce, in the limit where the spacetime approaches type D, to the Kinnersly tetrad 
used in Jj] for the BH background. In this tetrad the Weyl scalar ¥4 carries information 
on outgoing gravitational radiation, and on ingoing radiation. Work is in progress to 
identify one specific and physically significant tetrad from the equivalence class, appro- 
priate for a generic numerical relativity code using an arbitrary ADM slicing. However, 
the method of is already applicable when using a null slicing, in particular to the 
Bondi metric used in Q], where non-linear oscillations of a BH were analysed. In the 
case of the Bondi metric the gravitational wave signal can be extracted using the news 
function y v (v is retarded time), directly related to the outgoing energy. In the linear 
regime, one expects ¥4 = — y )VV . Thus this case is ideally suited to test our method, 
since we can compare |6] the news function obtained directly from the Bondi metric 
with that obtained via ¥4 in the Quasi-Kinnersley frame. Fig. El shows the comparison 
of ^4 with 7 VV , where in our axisymmetric case 2y = h^g, i.e. the other polarization 
fiQ^ vanishes. Clearly the agreement is excellent at late times, as predicted, with an error 
A = ^4 + y vv \ ~ 1(T 6 at v 4 = 80. 




CONCLUSIONS 

Here we have shown how the method obtained in fl, Eh works, reproducing very well 
results obtained in ED, when applied in the context of a code using the characteristic 
approach @]. Work is in progress to find a method to choose one particular Quasi- 
Kinnersly tetrad out of the general class in order to properly extract gravitational wave 
templates from any numerical relativity code using an arbitrary ADM slicing. 

For a rotating NS, we have shown that the deviation of the spacetime from Petrov 
type D is always very small, with 1 — \S\ rapidly decreasing out of the equatorial plan 
and with increasing distance from the star. 
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